Spatial separation of large dynamical blue shift and harmonic generation 
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We study the temporal and spatial dynamics of the large amplitude and frequency modulation 
that can be induced in an intense, few cycle laser pulse as it propagates through a rapidly ionizing 
gas. Our calculations include both single atom and macroscopic interactions between the non-linear 
medium and the laser field. We analyze the harmonic generation by such pulses and show that it 
is spatially separated from the ionization dynamics which produce a large dynamical blue shift of 
the laser pulse. This means that small changes in the initial laser focusing conditions can lead to 
large differences in the laser frequency modulation, even though the generated harmonic spectrum 
remains essentially unchanged. 

PACS numbers: 



In their pioneering experiments reported in 
Hentschel et al demonstrated the generation of single 
attosecond XUV pulses by 7 fs, 750 nm laser pulses inter- 
acting with a neon gas. As a first application, these XUV 
pulses were used to probe the sub-cycle time dependence 
of the laser electric field emerging from the neon gas. A 
surprisingly large dynamical blue shift was observed in 
this way, with a maximum value of approximately 35% 
of the laser frequency, a result which was not well repro- 
duced by theory. 

Ionization is intrinsic to all highly non-linear laser mat- 
ter interactions. It contributes free electrons to the in- 
teraction medium which induces a temporal and spatial 
variation in the refractive index during propagation. This 
causes both self-phase modulation and defocusing of the 
field 0. The experimental results in 1] raise a number of 
questions: i) Can such a large frequency modulation be 
caused by ionization driven self-phase modulation? ii) If 
so, how is it consistent with the observed harmonic cutoff 
energy around 90 eV, which indicates a relatively mod- 
erate interaction intensity around 5 x 10^^ W/cm^ f^? 
And iii) what is the role of the ionization dynamics in 
creating conditions favorable for the formation of single 
attosecond pulses? 

In this paper we answer these questions by analyzing 
the temporal and spatial ionization dynamics of an in- 
tense few-cycle pulse as it propagates through a neon gas 
cell and generates harmonics Q, using parameters simi- 
lar to those of the experiment in 1^ . We show that in the 
beginning of the gas cell a large laser frequency modula- 
tion builds up rapidly when the intensity is still high, and 
then slows down toward the end of the cell because the 
intensity is reduced by defocusing. We find that, quite 
generally, the harmonic generation is spatially separated 
from the frequency modulation of the driving field since 
the XUV radiation is built up predominantly in the sec- 



ond half of the medium. Furthermore, we show that the 
frequency modulation is very sensitive to small changes in 
the initial focusing conditions of the laser beam, whereas 
the intensity of the beam after propagation and the har- 
monic spectrum is not. Therefore, even though the XUV 
pulse is an ideal probe of the laser field after its interac- 
tion with the gas, the harmonic spectrum and in particu- 
lar its cutoff energy are poor indicators of the magnitude 
of the ionization driven dynamics. Finally we show that 
the ionization dynamics lead to a driving pulse at the end 
of the medium which has been spatially and temporally 
reshaped in a way that facilitates the generation of single 
attosecond XUV bursts. 

Our description of the intense-laser matter interaction 
includes both the response of a single atom to the laser 
pulse, and the collective response of the macroscopic gas 
to the focused laser beam. We solve the wave equation for 
the driving and the generated fields by space marching 
through the gas, using the slowly evolving wave approx- 
imation (SEWA) At each plane in the propagation 
direction we calculate the single atom dynamics by 
solving the time dependent Schrodinger equation. The 
non-linear atomic response is then used as the source 
terms in the wave equation for marching to the next plane 
in z. This approach is described in detail in 

For the propagation of the linearly polarized driving 
field Ei{t) we include two source terms due to the ion- 
ization of the medium 0: 
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where e and mg are the electron charge and mass, Ne{t) is 
the electron density, j{t) is the ionization rate, and Ip is 
the atomic ionization potential. All the time-dependent 
quantities are also functions of the cylindrical coordinates 
r and z. The plasma oscillation term Jp{t) gives rise to 
a spatial and temporal variation of the refractive index 
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FIG. 1: (Color online) (a) Time dependent ionization proba- 
bility of a neon atom exposed to a 750 nm, 7 fs laser pulse with 
a peak intensity of 10^^ W/cm^. Pvoi{t) is shown with solid 
line, Padk is shown with dashed line, (b) Intensity depen- 
dent correction factor for three different pulses: 7 fs duration, 
cosine carrier (open circles), 5 fs duration, cosine carrier (open 
triangles), 5 fs duration, sine carrier (filled circles). 



which causes defocusing and self-phase modulation. The 
absorption term Jabs {t) describes the loss of energy from 
the laser field due to the ionization of the medium. This 
term is small for all the cases discussed in this paper. The 
source term for the harmonic radiation is given by the 
non-linear polarization field, which is proportional to the 
single atom time-dependent dipole moment, calculated 
using the strong field approximation and the density 
of neutral atoms. 

To describe the short pulse ionization dynamics cor- 
rectly it is crucial to accurately calculate Ne{t) and ^{t) 
with sub-cycle precision. Our calculation of Ne{t) orig- 
inates in a numerical solution of the TDSE within the 
single active electron approximation l6|. We define the 
ionization probability Pyoi{t) from the probability den- 
sity of the wave function outside of a small volume around 
the ion core, which can be continuously evaluated during 
the calculation. In Fig. ^a) we show Pvoiif) (solid line) 
for a 750 nm, 7 fs driving pulse with a peak intensity 
of 10^^ W/cm^. Because of its short duration, the ion- 
ization probability at the end of the pulse is only about 
15%. 

Ideally, we would directly couple our numerical solu- 
tion of the TDSE to the solution of the wave equation 
and use Pvoi{t) to describe the time dependent ioniza- 
tion. Currently, we do not have the computational ap- 
paratus to do so. We also cannot directly use instan- 
taneous tunnel ionization rates proposed by Ammosov, 
Delone, and Krainov (ADK) 7] since these give rise to 
much larger ionization probabilities at high intensities 
than the numerical TDSE solution. However, we find 
that for a given peak intensity, the ionization probability 
PADK{t) calculated from ADK rates differs from Pvoi{t) 
by only a constant factor /3, as long as the intensity is be- 
low the (ADK) saturation intensity. We determine P as 
the ratio between Pyoi{t) and PADxif) at the end of the 
laser pulse. The dashed line in Fig.^a) shows /3PADK{t) 
which is in excellent agreement with Pyoi{t). 

We next calculate P{Io) by comparing PADK{t) and 
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FIG. 2: (Color online) Spatiotemporal dynamics of frequency 
modulation of 750 nm, 7 fs pulse propagating through 3 mm 
of neon gas. See text for details. 



Pvoi{t) for many different peak intensities, /q. This func- 
tion is shown in Fig. Hjb), open circles. We have cal- 
culated /3{Io) for different pulse durations and find that 
it depends only weakly on the duration and the abso- 
lute phase of the driving pulse. Examples are shown 
in Fig. njb) for driving pulse durations of 5 fs and 
two different carrier envelope phases. Finally, we use 
P{Io)PADK{t) to calculate the source terms A^e(^) and 
j{t) for each point in the non-linear medium, where Iq 
is the peak intensity of the driving pulse at that point. 
The insensitivity of P{Io) to the duration and phase of 
the driving pulse ensures that this is justified even as the 
pulse changes shape and phase during propagation. 

Fig. 121 illustrates the spatial and temporal dynamics of 
the frequency modulation experienced by a laser pulse 
during propagation through a rapidly ionizing neon gas. 
As in |l|, the incoming laser pulse has a wavelength of 
750 nm and a duration of 7 fs. The laser beam has a con- 
focal parameter of 4.2 cm and its focus is in the center of 
a 3 mm long neon gas jet with a density of 5 x 10^^ cm~^. 
In the absence of the non-linear medium the peak inten- 
sity in the focus would be 9 x 10^^ W/cm^. 

In[21^a) we show the on-axis electric field of the laser 
pulse after propagation through the neon gas. Its time 
dependent frequency, shown in (b), is calculated from the 
separation between consecutive peaks and zeros of the 
electric field in (a). The pulse exhibits a large frequency 
chirp with a shape characteristic of self phase modula- 
tion. The maximum blue shift which occurs near the 
peak of the pulse is very large, approximately 22%, al- 
though still smaller than the 35% shift that was reported 
in 1]. This will be discussed in more detail below. 

Fig.[2fc) shows how this large blue shift is accumulated 
during propagation. The dashed line (left axis) shows 
the on-axis energy density and the solid line (right axis) 
shows the blue shift, as functions of the propagation dis- 
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tance z. At each point in the blue shift has been found 
as the maximum value of the instantaneous frequency as 
plotted in (b). At the beginning of the gas where the 
intensity is high, defocusing is very strong. The energy 
density is reduced by almost a factor of two by the end 
of the gas. This intensity variation in turn controls the 
build-up of the blue shift which accumulates rapidly in 
the first half of the gas and saturates towards the end. 

The solid line in Fig. Etd) shows the radial variation 
of the blue shift at the end of the gas (left axis). The 
blue shift is largest on axis and decreases as a function 
of r. The dashed and dot-dashed lines (right axis) show 
the radial variation of the laser intensity at the begin- 
ning and at the end of the medium, respectively. The 
radial variation of the blue shift closely follows that of 
the incoming laser field, in agreement with the result in 
(c) that the blue shift is predominantly generated at the 
beginning of the gas. Experimentally, the laser field is 
probed by the XUV pulse after refocusing both beams in 
to a second gas jet by a mirror placed several meters from 
the first jet 1]. The XUV beam in the second medium 
is much smaller than the laser beam and therefore only 
probes its on-axis frequency modulation P. In principle, 
the radial variation of the dynamical blue shift could be 
probed by slightly misaligning the two beams. 

The ionization driven dynamics of the laser pulse has 
important consequences for the harmonic generation, and 
therefore for the formation of attosecond XUV pulses as 
these are synthesized from a range of harmonics near the 
cutoff 1]. Fig. [2l shows the radially integrated harmonic 
spectrum at the end of the neon gas. As a result of the 
frequency shift of the laser pulse, the harmonic struc- 
tures in the spectrum are not at odd multiples of the 
incoming laser frequency. As a result of the defocusing 
of the laser beam, the harmonic spectrum exhibits two 
different cutoffs. The high energy cutoff around 170 eV is 
determined by the peak intensity of the incoming beam 
(?^ 9 X 10^^ W/cm^) P, whereas the dominant low en- 
ergy cutoff around 90 eV corresponds to the reduced peak 
intensity in the second half of the medium. In an experi- 
ment, it is likely that only the low energy cutoff would be 
observed as the high energy cutoff is orders of magnitude 
weaker [Tq| . 

The insets show how the XUV radiation around 90 eV 
(left inset) and above 155 eV (right inset) build up during 
propagation. The highest energies are only generated 
over a short distance in the beginning of the medium 
and are then reabsorbed through the remainder of the 
gas. In contrast, the 90 eV photons are generated all 
through the medium. The saturation of this signal at 
the end of the medium is due to phase matching. For 
longer propagation lengths one would observe a periodic 
increase and decrease of the yield. The total energy in a 
5 eV range around 90 eV is approximately 10 p J. 

Fig.|31shows that the harmonic spectrum and its cutoff 
energy are poor indicators of the propagation dynamics 
of the driving field. In a rapidly ionizing, long medium 
the harmonic generation is spatially separated from the 
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FIG. 3: (Color online) Radially integrated harmonic spectrum 
at end of the neon gas. The insets show the build-up along 
the propagation direction of the XUV radiation around 90 eV 
and 155 eV (left and right inset, respectively). 



defocusing and self-phase modulation experienced by the 
laser field. In particular, both the yield and the pho- 
ton energy of the spectral cutoff are essentially decoupled 
from the dynamical blue shift of the laser pulse. 

This decoupling is further demonstrated in Fig. ^ We 
show results of a calculation where the focusing con- 
ditions of the incoming laser beam have been slightly 
changed compared to Fig. [2l while its energy has been 
kept constant. The new beam has a tighter focus with 
a confocal parameter of 3 cm and a peak intensity of 
12.6 X 10^^ W/cm^. This corresponds to decreasing the 
beam waist by less than 20%. 

Fig. 2la) compares the propagation dynamics of the 
two laser beams, showing the z-dependence of the on-axis 
energy density and the on-axis maximum frequency 
shift, as in Fig. ISl^c). The higher initial intensity gives 
rise to more free electrons, which in turn cause stronger 
defocusing. This means that the two beams have almost 
identical intensities at the end of the medium, and 
therefore give rise to very similar harmonic spectra, as 
shown in Fig. ^b). However, the peak dynamical blue 
shift of the more tightly focused beam is much higher 
than before and reaches almost 35% by the end of the 
medium. This answers the first two questions posed 
in the introduction: i) Ionization induced self phase 
modulation can indeed induce a 35% frequency shift 
of the laser pulse in conditions very similar to those in 
1], and ii) because of defocusing, the harmonic cutoff 
energy is essentially decoupled from the magnitude of 
the blue shift. 

To answer question iii) about the influence of the ion- 
ization dynamics on the formation of attosecond pulses. 
Fig. ISIJa) shows the three dimensional spatiotemporal 
profile of the intensity of the laser pulse at the end of the 
neon gas. We have used laser parameters as in Fig. 13 and 
the absolute phase of the input electric field is given by a 
cosine driver. The strong temporal and spatial reshaping 
of the laser beam that takes place in the beginning of the 
medium results in a broad, divergent beam with a large 
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FIG. 4: (Color online) We compare the build-up of the fre- 
quency modulation and the harmonic spectrum generated by 
two laser beams with the same energy but slightly different 
focusing conditions. Thick lines show results of a more tightly 
focused beam than that used in Figs. 2 and 3 (thin lines). 



radial phase variation, in which the peak of the pulse oc- 
curs at different times for different radial positions. In 
particular, this means that the pulse on axis has been 
shortened to a duration of less than two optical cycles 
(5 fs) compared to the initial 7 fs duration. 
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FIG. 5: (Color online) Spatiotemporal profile of the laser 
pulse at the end of the neon gas. In (a) we have used condi- 
tions as in Fig. El and (b) show a laser pulse that gives rise 
to the same harmonic cutoff energy around 90 eV, but has a 
lower initial peak intensity (5 x 10^^ W/cm^). 

The XUV radiation near the cutoff energy, which is 
used to synthesize attosecond pulses, is only generated at 
the highest intensities. The laser pulse shown in Fig.[5Ja) 



will give rise to one XUV burst on axis 0, and one or 
several bursts off axis half a cycle later. In the nearfield, 
the XUV time profile will therefore in general consist of 
several attosecond bursts. However, the off-axis XUV 
radiation is generated by a strongly divergent wavefront 
and is much less collimated than the on-axis XUV radi- 
ation. A spatial filter placed in the farfield can therefore 
also act as a temporal filter, selecting only the single, col- 
limated attosecond XUV burst generated on or close to 
the axis [llil. 

To demonstrate that the ionization dynamics play 
a crucial role in producing the laser pulse shown in 
Fig. [SJ^a), we show in Fig. [5Jb) the nearfield profile of a 
laser pulse with a lower initial intensity (5 x 10^"^ W/cm^ 
and a confocal parameter of 6 cm). This pulse gives rise 
to the same spectral cutoff around 90 eV. The absence 
of spatiotemporal reshaping of this pulse makes it much 
less ideal for generation of single attosecond bursts: there 
is no shortening of the pulse on axis, and XUV bursts 
generated in consecutive half-cycles of the field will have 
similar divergence properties because of the smooth ra- 
dial structure of the beam. 

In summary, we have analyzed the propagation and 
harmonic generation dynamics of an intense few-cycle 
laser pulse. We showed that the frequency modulation 
induced by rapid ionization in the non-linear medium 
can be very large, and can indeed reach a value of 35% 
in conditions very similar to those in 1]. We discussed 
that whereas the frequency modulation is very sensitive 
to the laser focusing conditions, the harmonic spectrum 
is less so. Most importantly, we showed that the ioniza- 
tion driven spatiotemporal reshaping of the laser beam 
is very important in creating conditions favorable for the 
generation of isolated bursts of attosecond XUV pulses. 
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